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Abstract 

The sodalite network is the edge-skeleton of the uniform tiling in Eu¬ 
clidean 3-dimensional space by Archimedean tetrakaidecahedra (truncated 
octahedra). We develop explicit expressions for its height (minimum net¬ 
work path length from some fixed to given vertex) and coordination (con¬ 
tent of network sphere of given height) functions. The final discussion 
should to some extent assist in motivating and signposting our proof strat¬ 
egy, in the course of ruminating on its potential generalisation. 
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1 Geometry and Symmetry 

The sodalite network is the edge-skeleton of the uniform tiling in Euclidean 
3-dimensional space by Archimedean tetrakaidecahedra (truncated octahedra). 
It may be constructed by ‘bitruncating’ vertices of the standard cubical tiling, 
so that octahedra forming around cubical vertices collide and in turn become 
truncated. 

After scaling up by a factor 4, the canonical cell centred (around a cubical 
vertex) at origin (0,0,0) has 24 vertices 

{( 0 , 1 , 2 ), ( 0 , 2 , 1 ), ( 2 , 0 , 1 ), ( 2 , 1 , 0 ), ( 1 , 2 , 0 ), ( 1 , 0 , 2 ), 

( 0 ,- 1 , 2 ), ( 0 , 2 ,- 1 ), ( 2 , 0 ,- 1 ), ( 2 ,- 1 , 0 ), (- 1 , 2 , 0 ), (- 1 , 0 , 2 ), 

( 0 ,- 1 ,- 2 ), ( 0 , — 2 , — 1 ), (- 2 , 0 ,- 1 ), (- 2 ,- 1 , 0 ), (— 1 , — 2 , 0 ), (- 1 , 0 ,- 2 ), 
( 0 , 1 ,- 2 ), ( 0 ,- 2 , 1 ), (- 2 , 0 , 1 ), (- 2 , 1 , 0 ), ( 1 , — 2 , 0 ), ( 1 , 0 ,- 2 )}, ( 1 ) 

generated by the octahedral group 0 3 acting on the canonical vertex O = 
(0,1,2). A useful rule of thumb is that P = (x, y, z) represents a network ver¬ 
tex just when none of the 6 sums or differences of pairs of integer components 
vanishes: 

x±y, y±z, z±xj^0 (mod 4). (2) 
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Each vertex P has four network neighbours, at Euclidean separation ^/2; the 
neighbours of O = (0,1, 2) comprise 

{(0,2,1),(0,2,3),(1,0,2),(-1,0,2)}. (3) 

A sector denotes any of 48 isomorplis of the canonical fundamental region 

U = {{x,y,z) | x > y > z > 0} (4) 

under the action of O 3. To any vertex P = ( x,y,z ) corresponds an isomorph 
P' = (x\ y', z') G U, where x\ y ', z' denotes the same bag of integers as x, y, z, 
apart from dropping signs and sorting. Remark that O = (0,1,2) / (0,0,0), 
the origin; also O U ! 

With translation symmetry 

T 3 : P —^ P T (2,2,2) of path length 3, associate its quadrant (union of 12 
contiguous sectors) 

{(cc, y, z) | x + y>0&y + z>0&;z + x> 0}; (5) 

T 4 : P —> P + (4, 0, 0) of length 4, its 8 -sector quadrant (sextant?) 

{{x,y,z) | x > \y\ & x > \z\}] (6) 

T 6 : P —> P + (4,4, 0) of length 6 , its 4-sector quadrant 

{{x,y,z) | y > \z\ & x > y + \z\}. (7) 

T 3 carries the canonical cell to one centred around a cubical centre; together 
with O 3 it generates the network symmetry group. 


2 Height Function 

Denote by ht(P) the height of vertex P , the minimum path length from O to 
P along network edges. 

We introduce a mildly indigestible function h(P) on vertices (revealed to 
equal height by Theorem in terms of h'(P') (equal to mean height over an 
O 3 orbit): 

Definition 1. For P' = (x', y’, z r ) G U, let 

h'{P')=x' +y'/2-s{P')/2, ( 8 ) 


where 


s(P') 


= (0,1, 0,-1) if y' mod 4 = 0,1,2,3, 

x (1,0,—1,0) if x' mod 4 = 0,1,2, 3. 


Now extending this to all sectors, 
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Definition 2. 


h(P) = h'(P') + 
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The sector offsets involved above are illustrated in Figure [TJ in the form of a 
section across the x axis, y running top to bottom, z left to right, and origin 
(beneath) centre. Remark how sectors with a common coordinate plane share 
equal offsets; also ( x,y,z ) —> (— x,y,z ) is a symmetry, (x,y,z) —>■ (x,z,y) an 
antisymmetry. 



Figure 1: Height offsets by sector: x out, y down, z right, origin central. 

Within its quadrant, a translation ‘respects’ the pretender height function: 

Lemma 3. 


h(T 3 P) = h(P) + 3 
h(T 4 P) = h(P) +4 
h(T e P) = h(P)+ 6 


for P € T 3 quadrant; 

for P e T 4 quadrant; 

for P £ Tg quadrant. 
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Proof. Within the T 3 quadrant, suppose without loss of generality that |x| > 
\y\ > \z\; then x > y > 0, and via Definition [2] 

/i(T 3 P) = {pc + 2) + (y + 2)/2 + s(x + 2, y + 2, z + 2)/2 
= x + y /2 + s(x, y, z)/2 + 3 = h(P) + 3. 

Similarly for T 4 , supposing that \y\ > \z\, 

h(T A P) = (a: + 4) + y/2 + s(x + 4, y, z)/2 

= x + y/2 + s(x, y, z)/2 + 4 = h(P) + 4. 

For T 6 , 

h(T 6 P) = (x + 4) + (y + 4)/2 + s(x + 4, y + 4, z)/2 
= x + y/2 + s(x, y, z)/2 + 6 = h(P) + 6 . 


□ 

Given a vertex pretending to sufficient height, a useable neighbourhood of 
it lies entirely within some quadrant: 

Lemma 4. If vertex P £ U with h(P) > 15, t/ien P together with its neighbours, 
and their images under the associated inverse translation T^ 1 , lie within the 
same quadrant. 

Proof. The components of any neighbour Q vary from those of P by — 1,0, +1. 
Referring to Equation [5] etc. — 

Case y + z > 6 : then 3 < y < x. At worst, 

P = (x, y, z), Q = (x, y - 1, 2 - 1), 

T 3 - 1 P= {x-2,y-2,z-2), T^Q = (x — 2, y — 3 ,2 — 3); 

all remain within the quadrant, including T 3 ^ 1 Q, since 

x + 2 / — 5>0&y + 2 — 6 > 0 & 2 : + x — 5>0. 

Case y + z < 6 : then 0 < y, z < 6 , and x > 12 via Definition [ 2 ] At worst, 

P = (x, y, z), Q=(x-l,y+ 1, z), 

T 4 - 1 P = (x - 4, y, z), T 4 - 1 Q = (x-5,y + l,z); 

all remain within the quadrant, including T 4 ^ 1 Q, since 

x — \y\ — 6 > 0 & x — |z| — 5 > 0. 


□ 


Theorem 5. ht(P) = h{P) for every network vertex P. 
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Proof. For ht(P) < 15 the assertion is verified via inspection of an inconve¬ 
niently extensive table. For ht(P) > 15 via a somewhat delicate induction: 
assume the result for all vertices R with ht(P) < ht(P). 

First suppose P £ U. P has some neighbour Q nearer to O , so that ht(Q) = 
ht(P) — 1; via Lemma[I]all P,T~ 1 (P),Q 1 T~ 1 (Q) lie in one quadrant, and via 
Lemma [3] and network symmetry 

ht(P) = ht(Q) + 1 = h(Q) + 1 = h(T~ x Q) + t + 1 = h(T~ 1 P) + t = h(P), 
where t = 3,4 for T = T 3 , T 4 resp. 

For any vertex P now, the procedure above may be applied to its sector, 
employing an appropriate symmetry from O 3 , and corresponding height offset 
from Definition [2] □ 


3 Spherical Content Function 

Consider now the content of a ‘sphere’ in the corresponding metric, also the 
‘vertex coordination sequence’ of the tiling: that is, the number of vertices S(n) 
at given height n from O. 

We introduce functions (there are more to come!) for the content of a ‘sphere’ 
of given height n, and of the ‘armillary sphere’ where it meets the coordinate 
planes, and of their restrictions to the canonical sector: 

Definition 6. 

S(n) = #(P | ht(P) = n); 

S(n) = #(P | ht(P) = n & xyz = 0); 

S 2 (n) = #(P G U | ht(P) = n); 

Si(n) = #(P € U | ht(P) = n & xyz = 0); 


n = 
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S{n) = 

1 
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40 

44 

48 

Sin) = 

4 

10 

20 

34 

52 

74 

100 

130 

164 

202 

244 

290 

340 


Figure 2: Table of network sphere content 


Lemma 7. There are finitely many constants aj , bj,Cj such that for all n, 

S{n) = Y^ a jSi(n+ j), 

0 

s ( n ) = y. bjSijn + j) + CjS 2 (n + j). 

3 3 
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Proof. Firstly notice that via Equation[2]the only tiling vertices on the boundary 
of a sector are those interior to a facet on a coordinate plane, since the other 
two planes are diagonal. 

Now via Definition [9] and Theorem [5J for n > 0, 

S(n) = 6Si(n) + 3Si(n — 1) + 3Si(n + 1) + Si(n + 2) + Si(n — 2) 


+ Si(n + 3) + Si(n — 3) + 2Si(n + 1) + 2S\{n — 1) 

+ 2S\{n + 2) + 2Si(n — 2); 

S(n) = 2(6S 2 (n) + 3S 2 (n — 1) + 3S 2 (n + 1 ) + S 2 (n + 2) + S 2 (n — 2) 
+ S 2 (n + 3) + S 2 (n — 3) + 2S 2 (?r + 1) + 2 5 2 (?z — 1) 

+ 2 5 2 (n + 2) + 2 5 2 (n — 2)) — 5(n); 


the final term corrects for the boundary being counted double, and is then 


substituted via the first equation. 


□ 


Theorem 8. 




Proof. P = (a:, y, z) will be restricted implicitly to vertices at height n in sector 
U. Also assume n > 6 , avoiding special values for n < 0. Lemma [3] is employed 
without reference, noting IA is a subset of all relevant quadrants. A leaning 
tower of further subsidiary functions follows: 

Firstly via P —> T 4 ~ l P from Equation [5] etc., 

Si(n) = #(P | 0 = 2 < y < x) = Si(n - 4) + S 0 (n); 

where via P —> T 6 ^ 1 P, 

So(n) = #(P | x — 4 < y <x) = So(n — 6 ) = constant, 

depending on n mod 6 , n mod 4. Hence the Si (12* + j) are polynomials linear 
in i, depending only on j = n mod 12 . 

Similarly via P —> T 3 - 1 P, 

S 2 (n) = #(P | 0 < z < y < x) = S 2 (n — 3) + Si(n); 
where via P T 6 ^ 1 P, 

Si(n) = #(P | 0 < ^ < 2) = Si(n — 6 ) + S 0 (n); 

where via P —> T 4 ^ 1 P, 

So(n) = #(P \ 0<z<2Szz<y<4:) = So(n — 4) = constant, 
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depending on n mod 3, n mod 6, n mod 4. Hence the S2(12i+j) are polynomials 
quadratic in i. depending only on j = n mod 12. 

For via LemmaQ] both S(12i+j), S(12i+j) are also such sets of polynomials; 
and by inspection of Figure [2J each set reduces (mysteriously) to the single 
polynomial in n > 0 shown. □ 


4 Discussion 

Coordination sequences for networks associated with lattices in Euclidean d- 
space are established in [3] and [2J for numerous classical cases, using general 
group-theoretic methods. However, an explicit expression for the coordination 
sequence S(n) in Theorem [ 8 ] — case d = 3 of Equation (3.43) in [3] — was 
designated conjectural, and has since been dubbed by one author ‘really tricky’. 

The approach used above relies on two ideas. One is intuitively obvious: 
for all vertices P lying sufficiently far from the height zero vertex O in the 
direction of translation T, translation must respect height, in the sense that 
ht(T P) = ht(P) + t for constant t. As a result, the sphere of height n + t is 
the union of translations of overlapping segments of the sphere of height n. The 
difficulty comes in quantifying regions in which this respectful behaviour can be 
guaranteed, for some neighbourhood of P sufficient to facilitate induction: it is 
overcome by overlapping adjacent quadrants so far, that P cannot avoid lying 
well inside (at least) one. 

Given the explicit expression Definition [2] for height, it would be feasible 
(though tedious) to compute content S(n) as the number of solutions of the 
compound linear equation h(x,y,z ) = n. However, it is not actually necessary 
to know ht(P) in order to find 5(n); only Lemma [3l that outside some finite 
initial region, translation respects height. 

References to ‘inspection’ are in practice made tongue-in-cheek. Frequent 
resort to a computer is inevitable if the proof is to be fully checked, starting with 
tabulation of ht(P) < 20 via some tree-search algorithm, and implementation 
of h(P), and extending either to compute S(n) by enumeration. 

Various features of the proofs merit further discussion. 

(A) The coordinates employed above for sodalite are special to dimension d = 3, 
and do not generalise to other spaces. Conventional coordinates for d-space 
sodalite (aka ‘honeycomb of permutohedra’) are sketched in pQ. 

(B) The proof of Lemma [7] and Theorem [H] involved establishing existence of 
moderately complicated linear combinations and polynomial sets; however, their 
actual coefficients are ultimately irrelevant. 

(C) Indeed, if only we could be certain in advance that the the final result 
in Theorem [5] would be a quadratic polynomial (or more generally, ultimately 
satisfied some linear recurrence with constant coefficients and known order), 
only minor routine computation would subsequently be required. 

(D) Although in this instance height ht(P) is represented by expression h(P) 
for all ht(P) > 0, content S(n) takes special values unless n = ht(P) > 1. In 
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other situations it may be nontrivial to establish a lower bound on height above 
which the corresponding functions behave well. 

(E) Larger still is the bound - ht(P) > 15 in Lemma 0] - required to ensure 
that height is respected by inverse translation, below which we are obliged to 
verify a result via inspection. Unfortunately this bound depends in a compli¬ 
cated fashion on the geometric interaction between translations, and may well 
grow rapidly with dimension. 

(F) The quadrant associated with a translation T would in general be the union 
of those sectors having an edge along its axis. This rule is subverted by T 3 
in Equation [5J which spills over into a further 6 adjacent sectors. The effect 
is to simplify the proof of Lemma [2 which otherwise would otherwise involve 
(dimension) d = 3 translations including Tg. 

(G) The proof of Lemma [7] would be considerably complicated by any necessity 
to consider vertices lying on lower-dimensional boundary elements (here edges, 
corner) of a sector. 

(H) Assignment of coordinates labelling vertices is a delicate matter, governed 
by availability of convenient generators for the symmetry group. In particular, 
it may be inadvisable to locate vertex O at origin (0, 0,0), or indeed to assign a 
vertex to the origin at all: instead the frame should reflect the full point-group 
of the network. 

All in all there appear to remain considerable obstructions to recasting our 
approach as an abstract theorem applicable to a class of networks, analogous to 
Theorems 2.4 and 2.9 of j3). In the meantime, the method has been successfully 
applied to the snub-square and knight’s move networks. 
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